Resonance fluorescence of a trapped three-level atom 
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We investigate theoretically the spectrum of resonance fluorescence of a harmonically trapped atom, 
whose internal transitions are A-shaped and driven at two-photon resonance by a pair of lasers, 
which cool the center-of-mass motion. For this configuration, photons are scattered only due to 
the mechanical effects of the quantum interaction between light and atom. We study the spectrum 
of emission in the final stage of laser-cooling, when the atomic center-of-mass dynamics is quantum 
mechanical and the size of the wave packet is much smaller than the laser wavelength (Lamb-Dicke 
limit). We use the spectral decomposition of the Liouville operator of the master equation for the 
atomic density matrix and apply second order perturbation theory. We find that the spectrum of 
resonance fiuorescence is composed by two narrow sidebands — the Stokes and anti-Stokes compo- 
nents of the scattered light — while all other signals are in general orders of magnitude smaller. 
For very low temperatures, however, the Mollow-type inelastic component of the spectrum becomes 
visible. This exhibits novel features which allow further insight into the quantum dynamics of the 
system. We provide a physical model that interprets our results and discuss how one can recover 
temperature and cooling rate of the atom from the spectrum. The behaviour of the considered 
system is compared with the resonance fluorescence of a trapped atom whose internal transition 
consists of two-levels. 



I. INTRODUCTION 



Cold trapped atoms are an ideal system for investigat- 
ing the quantum properties of the mechanical effects of 
photon-atom interaction. Pioneering experiments with 
optical lattices and ion traps have allowed to measure and 
characterize several properties of the scattered radiation, 
thereby gaining insight into the dynamics of the driven 
atoms and, in particular, of the mechanical effects of light 
on the atomic center-of-mass motion Recently, m 

experiments with single trapped ions it has been possible 
to measure with high precision the elastic component of 
the light scattered by a driven dipole 3, 4] , and to observe 
and characterize the Stokes and anti-Stokes components 
due to the harmonic motion in the trap Q , thereby con- 
firming theoretical predictions |^ . Lately, the prop- 
erties and the manifestation of the mechanical effects in 
the light scattered by these systems is experiencing re- 
newed interest in several experiments, investigating the 
coupling of radiation with single atoms and ions in opti- 
cal resonators 
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In this work, we investigate the spectrum of resonance 
fluorescence of a harmonically trapped atom, whose in- 
ternal transition is A-shaped, and which is cooled by 
two lasers tuned at two-photon resonance. This configu- 
ration is peculiar, since photon emission arises only due 
to the mechanical effects in the photon-atom interaction. 
In fact, when the coupling between internal and external 
degrees of freedom can be neglected (e.g. for copropa- 
gating laser bear ns) . this system exhibits coherent pop- 
ulation trapping 0, 0, 0| : The electronic stationary 
state is a stable coherence, that docs not absorb pho- 
tons due to destructive interference between the dipole 
excitation paths, leading to no emission of photons at 
steady state. In contrast, for laser configurations where 
two photon processes are Doppler-sensitive, internal and 



external degrees of freedom are coupled. If the atomic 
center-of-mass motion is confined in a steep potential, 
such configuration may allow for laser-cooling to the po- 
tential ground state Here, we study the spectral 
properties of the radiation scattered by the atoms in the 
final stage of the laser-cooling dynamics. 

Our theoretical analysis considers the quantum dy- 
namics of the internal and external degrees of freedom 
of the atom. It is based on the perturbative expansion of 
the atomic dynamics in second order in the Lamb-Dicke 
parameter, i.e. in the ratio between the size of the wave 
packet over the laser wavelength . We extend previ- 
ous theoretical investigations [3, |^ , which analysed the 
Stokes and anti-Stokes components of the radiation scat- 
tered by a trapped dipole. In those studies these com- 
ponents dominate over the MoUow inelastic spectrum of 
the bare dipole '2l| , which is mainly due to photon scat- 
tering at zero order in the Lamb-Dicke expansion. In our 
case, the spectral component of the bare three-level atom 
disappears due to destructive quantum interference, and 
additional features emerge, which allow further insight 
into the coupled dynamics between the internal and ex- 
ternal degrees of freedom of the driven atom. We analyse 
each spectral component, and discuss how to extract in- 
formation from these results about the atomic dynamics 
at steady state. 

This work is organized as follows. In section II the 
system is described and the dynamics is discussed qual- 
itatively. In section III we present the theoretical de- 
scription and evaluate the spectrum using the formal- 
ism of 1^, Hf. We apply perturbation theory combined 
with th e sp ectral decomposition of the Liouville opera- 
tor [2^l23l |. and calculate contributions to the spectrum 
at higher orders in the Lamb-Dicke expansion. In sec- 
tion IV the results are summarized and compared with 
the picture for the dynamics presented in section II. Fi- 
nally, we discuss our results in comparison with the spec- 



trum of a two-level transition driven by a plane and by a 
standing wave, as evaluated in 0,0. In the Appendices, 
several details of the theoretical derivation are reported. 
The reader who is not interested in the theoretical details 
can discard section III without loss of coherence in the 
presentation. 



II. MODEL AND QUALITATIVE 
DESCRIPTION OF THE DYNAMICS 



a) 




ki ■ cos (/) 



k2 ■ cos (1)2 



We investigate theoretically the spectrum of resonance 
fluorescence of a driven three-level atom, whose center- 
of-mass motion is confined by a harmonic potential, as 
depicted in Fig. 1. For simplicity, we consider one- 
dimensional motion along the a;-axis. The relevant elec- 
tronic transitions are arranged in a A-configuration, com- 
posed of two stable or metastable states 1 1 } and | 2 ) and 
an excited state | 3 ) . The transitions \j) | 3 ) are 
dipoles with moments d^- and linewidth 7j {j = 1,2), such 
that the linewidth of the excited state | 3 ) is 7 = 71 -I- 72 . 
The atom is driven by the bichromatic field E — Ei -f E2, 
with 



b) 



c.c. 



(1) 



where x denotes the atomic center-of-mass position. 
Here £j and are amplitude and polarization of the field 
modes at the optical frequency wl, j with wave vector kj , 
and (pj is the angle between the laser wave vector and the 
X-axis. The components Ej drive the dipoles dj and are 
tuned by the same detuning S from resonance, such that 
the states 1 1 } and | 2 ) are resonantly coupled by two pho- 
ton processes. A detector monitors the light scattered at 
the angle ip with respect to the x-axis, thereby measuring 
the spectrum of the intensity. 

Throughout this work, we investigate the manifesta- 
tion of the mechanical effects of the interaction between 
light and atom in the spectral signal. The system is in the 
regime where the size of the center-of-mass wave packet 
Aa; is much smaller than the wavelength of the incident 
radiation (Lamb-Dicke regime), and the laser field E 
cools the motion ,183 1^ . 

In the Lamb-Dicke regime, the dynamics of the driven 
atom can be described by a hierarchy of processes at the 
different orders in the ratio Ax/Xl, which accounts for 
the effects of the field spatial gradient over the center- 
of-mass wave packet. At zero order in Ax/Xl, internal 
and external degrees of freedom are decoupled, and the 
internal stationary state of the atom is the dark state 



(2) 



with the Rabi frequency fij — dj ■ ejEj/h, which we as- 
sume to be real. In this limit, at steady state the density 
matrix of the atom is the product pDfi of the density 
matrix for the external degrees of freedom fi and for the 
internal degrees of freedom pn = |V'd)(''/'d|- 




FIG. 1: a) Geometry of the lasers and position of the detec- 
tor with respect to the axis x of the atomic center-of-mass 
motion. We denote with ki cos (pi and k2 cos 02 the projec- 
tions of the laser wave vectors on the a;-axis. The detector D 
records the light scattered at an angle ^p with respect to the 
motional axis, b) Relevant electronic transitions. The stable 
or metastable states | 1 ) and | 2 ) are coupled by dipole transi- 
tions to the excited state | 3 ) . The lasers drive the transition 
I J ) ~* I 3 ) with Rabi frequency ilj and both are tuned from 
resonance by 5. The state | 3 ) decays with rate 7j into | j ) 
= 1,2). 



At first order in Ax/Xl the state IV'd) becomes unsta- 
ble due to the spatial gradient of the field over the finite 
size of the wave packet: Therefore, at steady state the 
density matrix of the atom is pst = Pd/^ + 0{Ax/Xl), 
where the correction 0{Ax/Xl) accounts for the pro- 
cesses due to the mechanical effects of the coupling be- 
tween light and atom. 

The dynamics of this system has been investigated 
in 01 in the context of laser-cooling. There, it has been 
characterized by two main time scales: A faster time scale 
To, for the scattering processes at zero order in Ax/Xl, 
where internal and external degrees of freedom are decou- 
pled, and a slower time scale Ti ^ Tg, where the effects 
due to the field gradient along the center-of-mass wave 
packet manifest. On the time scale Tq the internal dy- 
namics accesses the dark state | ipD ) , and the atom ceases 
to scatter photons. On the time scale Ti, the atom ab- 
sorbs light that is out of phase with the laser field due 
to the harmonic motion, thereby leaving the dark state 
and undergoing transitions that change the vibrational 
state. Then, hght is scattered at zero order in Ax/Xl 
and the atom reaccesses the dark state. In other words, 
on a coarse-grained dynamics the internal state of the 
atom is pn, whereas the correction term 0{Ax/Xl) in 
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Pst accounts for the processes which couple internal and 
external degrees of freedom, and give rise to the spectrum 
of emission. 

In Fig. 2 two spectra of emission of the dipole di |23| 
are shown for (a) Ax/Al « 7 x lO""* and (b) Ax/Xl ~ 
4 X 10"'^. The most striking feature is the visibility of 
the sidebands of the elastic peak, namely the two signals 
at ujl,i ± i^, compared to the rest of the spectrum. These 
motional sidebands are Lorentz curves of equal height, 
and their functional dependence on the frequency is plot- 
ted in the insets of Fig. 2. They originate from Raman 
scattering processes, where the initial and final state is 
I'j/'d) and the vibrational state is changed by one phonon. 
These processes occur on the time scale Ti, which also 
determines the linewidth of the resonances. 

The broad signals of the spectrum correspond to the 
MoUow — type inelastic spectrum [2^ |2^ and their typi- 
cal height is orders of magnitude smaller than the height 
of the motional sidebands for Ax/Xl ~ 7 x 10"**, while it 
is comparable for Ax/X w 4 x lO^''. We remark that the 
intensity of the elastic peak is at higher order in Ax/X^. 
In the following, we evaluate and discuss the spectrum in 
detail. 
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FIG. 2: Spectrum of resonance fluorescence S{lj) (in arbitrary 
units) as a function of the frequency to — ujl,i in units of i/. 
The parameters are ^1=^2= 8.5v, 7 = lOu, 71 =72, 
S — 35^, (j>\ = 0, <j)2 ~ 7!", corresponding to (n) — 0.005. In 
a) rji — ri2 — 0.01, in b) 7?i = 772 = 0.05. In both figures, the 
insets show the details of the curves at oj = u!l,i ± v. The 
contributions to the spectrum are plotted separately: The 
solid line corresponds to the signal of the Stokes and anti- 
Stokes components, the dashed line to the signal due to the 
MoUow-type inelastic spectrum. 



III. EVALUATION OF THE SPECTRUM OF 
RESONANCE FLUORESCENCE 

Let the detector measure the radiation scattered by the 
dipole di |23| . In the far-field the spectrum at frequency 
Lu is determined by iS(w) = X'S'(w), where x collects all 
prefactors which do not depend on uj. All results are 
rescaled by this common factor, which contains the dipole 
radiation pattern and is therefore a function of ip. At a 
different angle than the laser propagation directions, the 
expression 

poo 

S{uj)=Re dTe-'^'^-'^'-''>^{D''{x,T)D{x,0)) (3) 
Jq 

contains the frequency dependence, where D{x, t) is the 
generalized dipole lowering operator for the transition 
1 ) ^ I 3 ) in the reference frame rotating with the laser 
frequency ujl,i- By means of the quantum regression the- 
orem, the two-time correlation function in Eq. is de- 
termined by the Liouvillian C defined in the master equa- 
tion dp/dt = Cp for the atomic dynamics [l^. Thus, the 
dipole lowering operator at time t in Eq. ^ is given by 
D{x,t) = D{x)e^* with D{x) = e-'^''"'°^'^\l}{3\. Here, 
the centcr-of-mass position x is an operator acting on the 
atomic external degrees of freedom. The average (•) is 
taken over the atomic density matrix pst at steady state, 
given by Cpst = 0. 

We evaluate S{llj), Eq. ©, by applying the spectral de- 
composition of the Liouville operator £ [22,12^ according 
to the secular equations 

- Xp\ 
P^C = Xp\ 

with eigenvalues A and right and left eigenelements p^ 
and /5^, respectively. The orthogonality and complete- 
ness of the eigenelements is defined with respect to the 
trace, such that Tr{/0'^ p"^} = Ja'.a, where any density op- 
erator p can be decomposed as p = P^Trjp^p} |28l |. 
We define the projectors onto the eigenspace correspond- 
ing to the eigenvalue A as — p^ ® p^ leading to 

cv^ = v^c = xr^. (4) 

Their action on an operator X is defined as V^X — 
p^Trjp^X}. By applying this formalism, we rewrite 
Eq. |(2J) as 

=ReV- — -Tt{D\x)V^D{x)p,,]. 

^ i{uj - ujl,i) - X 

(5) 

Hence, the spectrum of resonance fluorescence is the 
sum of Lorentz and/or Fano-like profiles, centered at the 
imaginary part of the eigenvalues of £, with width given 
by the real part of A. The exact evaluation of the spec- 
trum for this kind of problem, with an infinite number of 
degrees of freedom, is a hard task. Nevertheless, an an- 
alytic solution can be found in the Lamb-Dicke regime. 
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In this limit, we make perturbation theory in the param- 
eter Ax/Ai on a spectral decomposition of of which 
the spectrum {A} and the respective eigenelements at 
zero order are known. 



A. Theoretical description 

The master equation for the atomic dynamics reads 

= Cp{t) = ^[H.p{t)] + K.p{t), (6) 

where H is the Hamilton operator for the coherent dy- 
namics and /C is the Liouvillian describing spontaneous 
emission. We decompose the Hamilton operator as 

H = H^,, + Ho + V{x), 

where Hq — ^SJ^j^i 2\^)i^\ gi'^^s the eigenenergies of 
the electronic states in the reference frame of the laser, 
with S = ojl,! — ll>i — ujl,2 — ^2, where ojj denotes the 
frequency of the transition | j ) — + | 3 ) . The term iJmcc 
describes the center of mass motion of the atom with 
mass M in a harmonic potential of frequency i^, 

o2 



1 



1 



(7) 



where x and p are the canonical conjugate variables de- 
scribing position and momentum of the atom, whereas 
a and are the annihilation, creation operators of 
a quantum of energy fiv^ respectively, such that x = 
.Jh/2Mv{a + at), and p = i.JhMv/2{a'^ - a). We 
denote with \n) the eigenelements of i?mcc, fulfilling 
Hraec\n) = hv[n + l/2)|n) with n = 0, 1, 2, . . . 

The term V{x) describes the coherent interaction of 
the atom with the lasers at the position x of the center 
of mass, 



h.c. 



(8) 



and the exponentials account for the recoil momenta 
:khkj COS 4>j of the atom when absorbing or emitting a 
photon. The operator /C in Eq. lO describes the sponta- 
neous decay, according to 

M0= -J(|3)(3|p(t)+p(i)|3)(3|) 

+ E ^^■|^')(3|P,W|3)(j-|, (9) 

where 71 + 72 = 7. Here, we have introduced 
1 

Pj{t)^ y"dcos6iA/'(cos6')e''=^""™"VW e~^'=^"'™"^ (10) 
-1 

which describes the momentum transfer Tikj cos 9 due to 
the photons spontaneously emitted at angle 9 with re- 
spect to the motional axis and with angular distribution 
M{cos9). 



B. Perturbative expansion in the Lamb— Dicke 
parameter 

In the Lamb-Dicke regime we can approximate 
exp (zbifcj-a; cos (/j) = (1 — ?y| cos^ (/5/2(2a'''a -I- 1)) ± 
irjj cos(/3(at -fa)-)- 0{rj^), where 



2Mv 



is the Lamb-Dicke parameter, corresponding to the ratio 
of the size of the oscillator ground state over the laser 
wavelength. Here, rjj is the parameter of the perturba- 
tive expansion, and it fulfills the relation rjj y^2{n) + 1 « 
Ax/Xl ^ 1. In second order in rjj expression Q has the 
form 



(11) 



where the subscript a = 0, 1, 2 indicates the correspond- 
ing order in the perturbative expansion. In order to eval- 
uate S'ct(w), we expand the operators C and D in power 
of rjj , yielding 



Do = 
Di - 

D2 = 



|1)(3|, 

— ifcia;cos?/;|l)(3|, 
~\klx^ cos^ V'|l)(3|, 



and 



Cip 
C2P 



[^^moc, P] ~ 



1 

i^ 

iE^V2x\p]+IC2p. 



-[Ho + ViO),p]+ICop 



(12) 
(13) 
(14) 



In Eq. (|12|l we have introduced the Liouville operators 
Ce and £/, which account for the external and inter- 
nal degrees of freedom. They are defined as CeP = 
1/ih [H^,,, p] and £ip = l/i?i [Ho + V{0), p] + /Cop. 

By determining Ca , we have used the expansion of the 
interaction term, Eq. V{x) = F(0) -f T^a; -h ^22:^/2 + 
0(ry3) with 



9" 

dx' 



x=0 



, a- 1,2 



(15) 



and the expansion of Eq. ^ , yielding 

7j 



ICoP^ f(2|j)(3H3)(j|-|3)(3|p-p|3)(3|) 



/C2P = /3 E ^J-^l \j){3\{2xpx~x^p~px^)\3){j\. 



Here, /? — J^^ dcos6' A/'(cos0) cos^0 is a constant. We 
remark that the first order term /Ci vanishes after aver- 
aging over the angles of emission 9. 
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From Eq. ^ it can be seen that internal and external 
degrees of freedom are decoupled at zero order. Hence, 
the eigenvalues of Co are 



^0 — + A/, 



(16) 



with Xe and A/ being the eigenvalues of Ce and Cj, 
respectively. The projector Vq in the corresponding 
eigenspace factorizes into the projectors Vj' and 'P^^ 
assigned to the internal and external degrees of freedom, 
according to 



' I ' E 



(17) 



The spectrum of £/ characterizes the dynamics of the 
three-level transition and the spectral properties of the 
radiation emitted by the bare atom. The eigenvalues of 
Ce take on the values Xe = i^J^, with £ = 0, ±1, ±2, . . . 
Each eigenspace at A^; is infinitely degenerate, and the 
corresponding left and right eigenelements are, for in- 
stance, jll^ — \n-\- t){n\, /^^ = + ^1- These eigenele- 
ments constitute a complete and orthonormal basis over 
the eigenspace at this eigenvalue. In particular, the pro- 
jector over the eigenspace at Xe = is defined on an 
operator X as 

^A.=i..^ = 5]^lTr^{AlX} (18) 

n 

= ^\n){n\X\n + e){n + i\, 

n 

where Tr e denotes the trace over the external degrees of 
freedom. 

At first order in the expansion in 77-,-, internal and ex- 
ternal degrees of freedom are coupled, and the degener- 
acy of the subspaces at eigenvalue A^; is lifted ■ 
The perturbative corrections to the eigenvalues Aq, to the 
eigenelements Pq, /5q, and to the projectors Vq are found 
by solving iteratively the secular equations at the same 
order p in the perturbative expansion, that is 



P P 



(19) 
(20) 



where p^ and are the a-order corrections to the 
eigenelements p^ and pg. The explicit forms up to sec- 
ond order are derived in Appendix A. In particular, 
Pq^° = pdP is the steady-state density matrix at zero 
order, where /x is the density matrix for the external de- 
grees of freedom in the final stage of the laser-cooling 
dynamics | ,3Ct |. and has the form 



1 



in) 



where 



(n) = Trja^a/i} 



(21) 



(22) 



is the average phonon number at steady state. 

By substituting the explicit form of the operators into 
Eq. we find that the zero- and first-order contribu- 
tions to the spectrum vanish (see discussion in Appendix 
B), yielding 



5(0;) = 52(0;) + 0(77^), 



with 



.9(A) 



(lu - ujl.i) - A' 



(23) 



(24) 



Here, g{X) is a complex-valued function, which we de- 
compose for convenience into g{X) — /(^^(A) + /'^^(A), 
with 

f^HX)=Tr{Dlv^Dopi}, 

f^'HX)^Tr{Dlr^Dop2}. 

Using p^, P2 and Vi, evaluated in Appendix A, and mak- 
ing use of relation H17|l . we separate the trace terms /'■^•' 
and Z^^-' into the product of the trace over the external 
and over the internal (Tr/{}) degrees of freedom. By 
applying the cyclic properties of the trace and the the 
completeness relation for the external degrees of freedom, 
J:^^ P^- = 1b, we find 

/^'H^) = ^EK'OTr^,{(P^'-[:r,p])a;} 

1 „ 1 



xTr, <j Dlrf' 



-Do 



Ao — ^E ^ A^ + Ci 



{{Vfpx)x] 



xTr, <j Dlr^' 



^1, 



-Do-r 



Ao — ^'e ^ A^ + £/ 
+Sy^,o^^rE{{V^^[x,^i])x} 



VipB 



[Vi,Pb] 



xTr/ { Dl 



Ao — Cj 

+6y^,oTTE{{V^^px)x} 



Vi,Vt'Do 



-VipB 



xTi'r < D, 



° Xo — Ci 



Vi,r^'Do ^ 



Xe + C, 



■[VuPb] 



where we have used the relation Tte{Pe'^ X} 
Sxe,oT^'^'e{X}. Analogously, 

Z^'^A) = -1^Sx,^oJ2[^te{{V^^Px)x} 



xTti <j dIv^'DoC] 



TrE{ir;-[x,p])x} 
Tr/ {dIv^'DoCJ^ 



X'e + Ci 



[VuPB_ 



-ViPD 



X'e + Ci 

^-^5x,,oTrE{px^}TTi{DlP^'DoCj'[V2,PD]} 
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From the properties of the terms Tte{X} one can al- 
ready see that the eigenelements and eigenvalues oi Ce 
contributing to the spectrum at second order are at 
\e = 0, ±ii/. 

In summary, the first non- vanishing contribution to the 
spectrum of emission is in second order in the pertur- 
bative expansion. It can be decomposed into the sum 
of curves, centered at the imaginary part of the eigen- 
values A of the Liouville operator C, and weighted by 
the factor g{X). The eigenvalues are here determined up 
to the second order correction, A = Aq + A2 -I- 0{r]j), 
whereby Ai = 0, as shown in Appendix A and in |2fl| |. At 
zero order in the perturbative expansion Xq = Xj + Xe, 
where A/ is the eigenvalue of the Liouvillian of a bare A- 
transition, while the only relevant external eigenvalues 
are A^; = 0, ±izy. 

Below, we analyze the spectrum in detail. For later 
convenience, we rewrite S{uj) = S'oi(w)-t-S'M(w)-|-S'sB('^), 
where Sci{io) is the contribution of the elastic peak, at 
A = 0, the term 



ing to the projectors in Eq. H18|) . giving 



5- 



(25) 



denotes the contributions at A/ 7^ 0, which we refer to as 
the MoUow-type inelastic component |2^ . The term 



E 



5(A) 



i(t^ - cjl.i) - A 



(26) 



represents the contributions at A/ = 0, A^; = ±i^, that we 
identify with the sidebands of the elastic peak, or Stokes-, 
anti-Stokes components of the scattered radiation. 

We remark that 52(1^) does not depend on the posi- 
tion of the detector: the terms containing the perturba- 
tive corrections Di, D2 and their adjoints do not con- 
tribute to S{<jj) in second order in rjj. Furthermore, the 
Lamb-Dicke parameters 771 , 772 appear always in the form 
77iCOS<^i, 772COS(/)2, since /C2Po^" oc 1C2P'd = (see Ap- 
pendix B). Therefore, the mechanical effects in second 
order are solely determined by the scattering of laser 
photons and not by recoils due to spontaneously emit- 
ted photons. This behaviour is due to destructive quan- 
tum interference at zero order in the Lamb-Dicke ex- 
pansion, implying that light absorption is a first-order 
process |l&.i31j. 



1. The MoUow-type inelastic spectrum 

At zero order in the Lamb-Dicke parameter, the eigen- 
values A/ 7^ determine the position and the shape of the 
contributions to the MoUow-type inelastic spectrum. As 
A appears in the denominator of S'm(w), the second-order 
correction A2 can be neglected. 

The trace terms over the external degrees of freedom 
are conveniently evaluated using the basis set correspond- 



TrE{i'PE'^fJ'X)x} = xl [5xE,i,^i(n) + 1) + SxE-i^{n)] , 
TTE{x^fi} = xl[2{n) + l], (27) 



^ and 



with Xq — ^hjlMv. Thus, at second order in the 
Lamb-Dicke expansion the eigenelements of Le deter- 
mining the spectrum are in the eigenspaces correspond- 
ing to As =0, ±ii^. Since 5m (w) is linear in these terms, 
this component of the spectrum scales with 77^ cos 
is linear in the average phonon number (tt,). 

This spectral component is constituted by the con- 
tributions of the signals centered at Im{A/} and at 
Im{Aj} ± V. The latter originate from the term /^^•'(A) 
in g(X). We illustrate this behaviour in Fig. 3, where 
the MoUow-type inelastic spectrum is shown for different 
values of the detuning (5 (and correspondingly of the aver- 
age phonon number at steady state (n)). The frequencies 
Im{A/} are marked with crosses on the frequency axis, 
where the arrows indicate the center-frequencies Im{Aj} 
of which only the sidebands are visible. In order to high- 
light this splitting, we have taken borderline parameters, 
such that all peaks are clearly resolved. For realistic pa- 
rameters, the sidebands on the left side of the spectrum 
are visible, as shown in Fig. 2(b). 

The curves at Im{A/} can be reproduced by evaluat- 
ing the spectrum of emission of a bare three-level atom, 
whose ground state coherence has a finite lifetime p5ll26l| . 
Thus, they can be identified with the spectrum of the 
photons scattered at zero order in the Lamb-Dicke pa- 
rameter. The signals centered at Im{Aj}±i/ are peculiar. 
They correspond to processes where photon scattering is 
accompanied by a change in the vibrational state. Look- 
ing at the corresponding eigenelements, it follows that 
they stem from the inelastic processes which take the 
atom out of the dark state. We remark that in Fig. 2 
the left pole Im{Aj} falls at the same frequency as the 
left motional sideband: In fact, the parameters have been 
here chosen, so that the frequency of absorption along the 
cooling transition coincides with the narrow resonance 
characterizing the excitation spectrum [TsL [l9l l3^ . 



The sidebands of the elastic peak 

The spectral contributions at A/ — 0, A^; = al- 
low for a compact analytic form. Only the term /^^^(A) 
contributes to g{X) in S'sb(w)- After some algebraic ma- 
nipulation, we write 

g{XE) = TrE {{V^E^ia + a^)fi){a + at)} |/(Ab)P, 



where 



/(A^) = ^Trj{Dlj-^^[VuPD]}. (28) 
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The explicit form H28|l is found by applying the relation 

OG 

(A — Ci)^^ = J dte^'-^^'"'^* . Using the quantum regres- 



sion theorem we arrive at 



/(Ab) = -i 



(29) 



with ~ flf + and where we have introduced 
rj = Xq (fci cos 01 — k2 cos (p2 ) ■ The explicit form of Ssb {^) 
is determined after evaluating the second order correc- 
tions A2 to Ao = ±iiy. These are found by solving the 
eigenvalue equations (|Aip and ljA2|l at Aq = ztiiy. In the 
subspace at Xe = '^iv, A/ = 0, after tracing over the 
internal degrees of freedom, they read 

+ s{-v) [a)jl^^a - aa^fi^^] + H.c, (30) 

where /x^^ are the right eigenelements of Eq. H30|l at the 

eigenvalue A2. The left eigenelements jl fulfill the cor- 
responding equation for the action to the left l23l| . 
The coefhcient siy) is given by 



1 



2Mi> 



= V 



1^2(^2 _^4j,(i^/2- + (5)) ■ 



(31) 



We rewrite s(±i') — Re[s(±z^)] -I- i Im[s(±zy)], and define 
A± = 2Re{s{^v)}. Substituting into Eq. (|Sn)l. we obtain 
the more familiar form [20l 1231 



A2/i^= = -w[ata,/i^^] (32) 
+A_ [2ail^^a^ - a'^afi^^ - fl^'^a^a] 
+A+ [2at^^^a - aa^/l^^ - p.^^aa''] , 

with 9 — Im[5(j^)] + Im[S'(— i^)]. The corresponding 
eigenvalues are solutions of Eq. (|A8|) , and have the form 



X2{N,£) = -ieD -{2N+\e\) (A_ -A+), 



(33) 



where the index = 0, 1, 2, . . . accounts for the removed 
degeneracy inside the eigenspace. The explicit form of 

the corresponding left and right eigenelements fi ' , /i^'^ 
can be found in |22j| . They form a complete and or- 
thogonal set with respect to the trace over the exter- 
nal degrees of freedom, Trsj/i ' fl'^ '^'} = 5N,N'5t,t' and 



'^N ii^ ® /i ' " = Ie- In particular, the eigenelements 
eigenvalue A^; = iii', such that 



jl^'^, jl^'^ form a complete basis over the subspace at 



We remark that the density operator given in Eq. (|21|l is 
right eigenelement at A^ = £ = 0, that is /i = pP''^ . Using 
this basis for evaluating the trace terms over the external 
degrees of freedom, we get 



5sb(u;) = Re^ 



\I{^e)? 



i{uj — ujl,i — (.1^) — A^'^ 



X TiE{{a + at)^^-^}TrB{A'^''(a + a^)ii}. 

By using the explicit form of the eigenelements in [2^ we 
find that only the terms at A^ = 0, £ = ±1 contribute to 
the sum, giving 



'S'sb(w) 
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So, (34) 



where 75 = A_ — A+ and sn is the height at the center 
frequency and has the form j33| 



So 



47sJrj4(rj2 _ 4iy2) 



(35) 



Here, we have defined 73 = 7s Z??^, which is at zero order 
in the Lamb-Dicke expansion. 

From Eq. one sees that both sidebands of the elas- 
tic peak have the same form, independent of the angle V' 
of the detector with respect to the axis of the motion. 
In particular, they have the same Lorentzian shape, as 
shown in the insets of Fig. 2, and are centered at the fre- 
quency ±{v + v), where is a shift in second order in the 
Lamb-Dicke parameter. We illustrate this effect in Fig. 
4, where we have chosen suitable parameter to show this 
small shift most clearly. Here, v can be identified with 
the a. c. -Stark shift arising from off-resonant coupling to 
other dipole transitions at different vibrational numbers. 
The width 73 of the sidebands is at second order in the 
Lamb-Dicke parameter and corresponds to the cooling 
rate [l9l |. The height sq is in zero order in the pertur- 
bative expansion. With some algebraic manipulations, 
using (n) = A_|-/7s, it can be rewritten as 



So^^{n){l + {n)). 



(36) 



N 



Thus, the motional sidebands are well distinguished 
compared to the MoUow-type inelastic spectrum for 
if < {n). 



3. The elastic peak 

The contribution at the eigenvalue A = corresponds to 
the elastic peak, i.e. to the coherent part of the spectrum. 
In this system its appearance is due to the mechanical ef- 
fects of light: In fact, at zero order in the Lamb-Dicke ex- 
pansion there is no photon emission at steady state. We 
evaluate the radiation scattered at this frequency starting 
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FIG. 3: Spectrum S{uj) in arbitrary units as a function of 
<^ — uiL,i in units of u. The figures are at different values 
of the detuning 5, for Q,i = 0,2 = IQv, 7 = 5i/, 71 =72, 
(?!>i = 0, 02 = TT, "Ui = r)2 = 10~*. (a) 5 = 0.5z^, corresponding 
to (n) = 30. (b) (5 = 3.5;^, corresponding to (n) = 3.8. (c) 
5 = IQu, corresponding to (n) = 1. In all figures, the crosses 
on the frequency axis indicate the positions of the frequencies 
ImjA/}. The arrows indicate the frequencies Im{Aj} of which 
only the sidebands appear in the spectrum. 



from expression (jsj, 

= 7r<5(c^-c^i,i)|Tr{i?V,t}P, (37) 

which is the well-known form of the elastic peak contri- 
bution [3^ . The perturbative expansion of D and pst can 
now be applied for evaluating the average dipole moment 
TrjDpst}- We find the first non-vanishing contribution 
at 0(772), such that Tv{Dp^t} - TrjDjpa} + Tt{d\pi} + 
0{r]j). This signal is due to the lowest order correc- 
tions of the Debye-Waller factor, exp(— r^^ cos^ </'i/2), on 
the transitions \4>'D,n) —f |3,n) — > lipDjU), and to the 



pj cos ip on the transitions \ipr,,n) — > 
3,rt ± 1) \ip-D,n). Thus, the intensity of the radia- 
tion scattered at the elastic peak is at fourth order in the 
Lamb-Dicke parameter and it depends on the angle of 
emission. The finite life time of the dark state suggests 
also a broadened signal at this frequency. Our analysis 
shows that such contribution is of higher order in the 
perturbative expansion. 



IV. SUMMARY OF THE RESULTS AND 
DISCUSSION 

The spectrum of emission of a trapped ion, whose inter- 
nal degrees of freedom constitute a A-transition driven 
at two-photon resonance, is a remarkable manifestation 
of the mechanical effects of light: In fact, at steady state 
photons are emitted due to processes where the vibra- 
tional degrees of freedom are excited by absorption of a 
photon. We have evaluated the spectrum using pertur- 
bation theory in the Lamb-Dicke parameter. According 
to our results, the signal of the emission spectrum is in 
second order in the Lamb-Dicke expansion. We classify 
the spectral features into three main contributions, which 
we summarize below. 

At the laser frequency lo^^i the spectrum exhibits a 
(5-peaked signal, visible for instance in Fig. 2, that we 
identify with the elastic peak. This signal is at fourth 
order in the perturbative expansion. This order of mag- 
nitude is understood, as the dipole moment at steady 
state scales with r/'j ~ (Ax/Xl)'^. In fact, this signal 
is due to Rayleigh scattering processes where the initial 
and final state is the dark state |V'd) and the vibrational 
number n is conserved. Thus, this excitation originates 
from the lowest order mechanical corrections to the Rabi 
frequency and it depends on the angle of emission tp. 

At the frequencies lul^i^v one observes two narrow res- 
onances. These are the motional sidebands of the elastic 
peak, the Stokes and anti-Stokes components of the scat- 
tered light. They correspond to Raman scattering where 
the initial and final internal state is | ipu ) and the vi- 
brational number is changed by one phonon. The curves 
are Lorentz curves, whose dependence on the physical 
parameters is given by Eq. and plotted in the in- 
sets of Fig. 2. Their width is in second order in the 
Lamb-Dicke expansion and corresponds with the cooling 
rate . The height of the curves at the center frequency 
scales with the average phonon number {n) according to 
{ti){1 + (n)), so that the total intensity emitted at this 
frequency is proportional to rf'{n). Finally, the center 
frequencies of the sidebands are shifted from the zero- 
order center frequency by a contribution v at second or- 
der in the Lamb-Dicke expansion: This corresponds to 
the a. c. -Stark shift of the ground-state coherence due to 
the off-resonant coupling with the excited state. Figure 4 
shows the shift v for one sideband. 

The other spectral features, visible for instance in Fig. 
2(b), can be identified with the MoUow-type inelastic 
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FIG. 4: Spectrum of resonance fluorescence S{ij) (in arbitrary 
units) as a function of the frequency to — ujl,\ in units of 
v. The parameters are <^i = 0, <^2 = ti", rji = rj2 = 0.05, 
f^i = S72 = 8.5/^, 7 = lOz^, 7i = 72, (5 = 15:^, corresponding 
to (n) = 0.2. In the onset, the Stokes component is shown: 
The vertical dashed hue indicate the position of the frequency 
uo ~ = —I'- In the inset, the whole spectrum is shown. 
The signal of the inelastic part is here two orders of magnitude 
smaller than the sidebands. 

spectrum. These can be decomposed into the sum of 
Lorentz curves, whose height scales with and is Un- 
ear in (n), whereas the width is in zero order in the 
Lamb-Dicke expansion. Part of these features can be 
reproduced by evaluating the incoherent spectrum of a 
bare A-transition driven at two photon resonance, whose 
ground-state coherence has a finite decay time [2^ |2g . 
Hence, their origin can be explained with photon scatter- 
ing at zero order in the Lamb-Dicke parameter, occur- 
ring once the atom has left the dark state. Nevertheless, 
this part of the spectrum exhibits also peculiar features, 
which cannot be understood in these terms. These are 
in fact curves which characterize the excitation spectrum 
of the bare A-atom [S^j , and which here appear shifted 
by the frequency :kv from their center-frequency. They 
thus describe scattering processes where the vibrational 
number is changed by one phonon. For saturating driv- 
ing fields, they can be interpreted as Raman scattering 
processes, where the initial state is the ground-state co- 
herence IV'd) and the final state is another dressed state 
at a different vibrational number state. The linewidth 
of the emitted photon is then the linewidth of the cor- 
responding dressed state transition, while the center fre- 
quency is the corresponding a. c. -Stark shift. Here, the 
center frequencies are shifted by the trap frequency 
since the dark state is excited only by processes chang- 
ing the vibrational number. 

Remarkably, in second order in the Lamb-Dicke expan- 
sion the spectrum S'(w) does not depend on the position 
of the detector, apart for the dipole pattern of emission. 
This is another consequence of the fact that at zero order 
in the perturbative expansion the atom at steady state 
is de coup led from radiation because of quantum interfer- 
ence [31j . 

Using these results, one can characterize the steady 
state of the motion. For instance, the measurement of the 



linewidth of the motional sidebands gives the cooling rate 
of the process. The phonon number (n) at steady state 
can be measured through the ratio between the heights 
at the center-frequency of the motional sideband and of 
one peak of the incoherent spectrum. In this way, one 
gets a simple equation at second order in (n) whose coef- 
ficients are determined only by the laser parameters. We 
remark that for larger values of (n) the visibility of the 
motional sidebands over the MoUow-type inelastic spec- 
trum increases. This is illustrated in the inset of Fig. 4. 

In summary, the spectrum we have evaluated allows 
to gain insight into the quantum dynamics and steady 
state of the atom interacting with light. Our results 
are in agreement with the dynamical picture presented 
in section II. This picture is based on a clear separa- 
tion between the two time scales Tq, Ti, on which also 
the validity of the perturbative expansion lies. Analyti- 
cal estimates and numerical checks of the validity of this 
coarse-grained dynamics have been presented in |19| . 

It is interesting to compare these results with the 
features found in the emission spectrum of a trapped 
two-level atom. In a two-level transition driven by a 
plane wave, the incoherent spectrum has a contribu- 
tion at zero order in the Lamb-Dicke expansion, as at 
this order the internal steady state of the system is 
characterized by non- vanishing occupation of the excited 
state. The features due to the mechanical effects man- 
ifest here in the motional sidebands. These are nar- 
row resonances, whose width is the cooling rate. How- 
ever, at a fixed detection angle the curves are Fano- 
like profiles, whose relative height varies with -0 (while, 
once integrated over the solid angle of emission, have 
Lorentz shape and are equal) |g • This behaviour is 
an interference effect between Raman processes at sec- 
ond order in the Lamb-Dicke expansion 1^: Given 
e) ground and excited states of the dipole transition, 
the Raman processes \g,n) — > |e,n) — > \g,n± 1) and 
\g,n) — > |e,n ± 1) — > \g,n ± 1) are of the same order 
and lead to the emission of the photon. They therefore 
interfere, and their interference signal (the height of the 
sidebands) is modulated by the emission angle. 

This behaviour disappears when the dipole is at the 
node of a standing wave: Then, the carrier transition 
g, n) — > |e,n) is suppressed and at this order only the 
transitions \g.,n) |e,n± 1) — > |5,n± 1) occur. Hence, 
both motional sidebands are Lorentz curves of equal 
shape, independently of the emission angle (which just 
affects the total height of the signal, according to the 
dipole pattern of radiation). Thus, in this case transi- 
tions in zero order in the Lamb-Dicke expansion vanish 
because of the spatial mode structure, while light scat- 
tering occurs because of the spatial gradient of the field 
intensity over the center-of-mass wave packet. Remark- 
ably, in our case transitions in zero order in the Lamb- 
Dicke expansion are suppressed because of quantum in- 
terference between dipole excitation paths, and photon 
scattering occurs due to the phase gradient of the field 
over the center-of-mass wave packet. 
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V. CONCLUSIONS 

We have presented a theoretical study of the spectrum of 
fluorescence of a trapped atom whose internal degrees of 
freedom are driven in a A-configuration at two photon 
resonance. In this system, the atomic emission at steady 
state is only due to the mechanical effects of the atom- 
photon interaction. The spectrum has been evaluated at 
second order in the Lamb-Dicke expansion, i.e. in the 
expansion of the size of the atomic wave packet over the 
laser wavelength. We find that the spectrum is charac- 
terized by two narrow resonances corresponding to the 
motional sidebands, i.e. the Stokes and anti-Stokes com- 
ponents, and by a MoUow-type inelastic spectrum, while 
the elastic peak is at higher order. Through these prop- 
erties, several information about the quantum dynamics 
and steady state of the driven atom can be extracted, like 
the cooling rate and the temperature, and the contribu- 
tions of the individual scattering processes can be identi- 
fied. Furthermore, for relatively large temperatures the 
sidebands of the elastic peak may be orders of magnitude 
larger than any other spectral signal, and the spectrum 
can be said to be solely composed of these two frequen- 
cies. 

Our results provide an interesting insight into the 
underlying physics of the mechanical effects of light- 
atom interaction, and may contribute to on-going experi- 
ments investigating and engineering the coupling of single 
trapped atoms and ions with electromagnetic fields. 



Ao — 'Co 



(Ai-£i) + (A2-£2^ 



Analogously, one finds the perturbative corrections to the 
left eigenelements solving equations 1^01) at second 
order. This in turn allows to evaluate the perturbative 
corrections to the projectors Vq. Using = + 
Pi ® Pq, we obtain 0] 



(1 - V^)Vi 



Ao — -Co 



Ao — L-O 



(A5) 
(A6) 



The equations for the corrections Ai, A2 to Ao are found 
by multiplying Eqs. (|XT)l . (|X2|l by p^ on the left and 
taking the trace. The resulting equations are 



Ai = Tr{p^/:ip^} = 0, 

A2 =Tr{p^/:2P^}-Tr{p;^(Ai 



(A7) 

^i)p^}, (A8) 



Tr{p^/:2Pa + Tr{/5^ A 



-^0 ~ -Co 



where we have used Eq. (|A3p and relation p^Co — XqPq- 
From Eq. (|A7|) it is visible that Ai = for all eigenvalues 
A: In fact, the Liovillian Ci couples subspaces at differ- 
ent Xe, but it vanishes inside a subspace at fixed Ae, 
V^LiV^ = 0. 
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The term at zero order in the perturbative expansion 
^"(^) = E 7 ^-^^"^^{dIv^^DoPo} = 

^ - t^L.l) - Ao L J 



vanishes, since 



APPENDIX A: PERTURBATION THEORY 

The equations (fl^ to solve iteratively in the perturbative 
expansion are 

Cop^ + Cip^ = Xop^ + Xip^, (AI) 
Cop^ + Cip^ + £2^0 = ^P2 + AiPi + A2P^,(A2) 

where p^ satisfy Cop^ = AoPo- For Aq = 0, po = poM, 
with p given in Eq. I|21|l . Equation (|A1() gives 



(1 - V^)Pi 



1-Vl 
Ao — -Cq 



(Ai-/:i)p^, (A3) 



where Vq is the zero-order projector onto the subspace 
at eigenvalue A, Vq = p^® Pq. Inserting (jA3|) in HA2|I we 
obtain 



(l-^0>2 



i-^r 

Ao — Co 



(A4) 



DpQ = pqD^ = 0, 



(Bl) 



as there is no excited state occupation at steady state in 
zero order. Analogously, in first order it can be shown 
that 

= Ev ^-^^{"^^{dIv^Dqpq] 

+ Tr [d\V^DqPq^ + Tr {^dIv^D^pq'^ 

+ Ti[dIv^DqpiY) =0, 

where each of the first three terms are equal to zero be- 
cause of relation (|Bip . The last term is equal to zero 
because here the position operator x occurs linearly. 

The non-vanishing contributions to the spectrum at 
second order are shown in Eq. (|24|l . All other terms van- 
ish. For most of them, this can be demonstrated using 
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(jBl|l . We would like to emphasize the disappearance of 
the contributions 

TT{Dlr^DoPo} = Tv{Dlv^Dipo} = 0, 
Ti{DlV^DoPi} = Tv{Dlv^Dipi} = 0, 
Tr{DlP^DoPo} = Tt{dIv^D2Po} = 0, 



which, together with ICpo = 0, imply that the spectral 
signal does not depend on the angle of emission up to 
second order. 
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